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Abstract In the Earth-Moon system, low-energy orbits are transfer trajectories 
from the earth to a circumlunar orbit that require less propellant consumption 
when compared to the traditional methods. In this work we use a Monte Carlo 
approach to study a great number of such transfer orbits over a wide range of initial 
conditions. We make statistical and operational considerations on the resulting 
data, leading to the description of a reliable way of finding "optimal" mission 
orbits with the tools of multi-objective optimization. 
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1 Introduction 

When facing the problem of transferring a spacecraft from an orbit around the 
Earth to one around the Moon, the option usually taken is that of a two-impulse 
Hohmann transfer. The past two decades have seen the discovery and develop- 
ment of a new family of orbits that exploit the A''-body dynamics to reduce the 
Av budget of interplanetary missions. These "low-energy orbits" are rather new 
but have already been used several times to reach the Moon and the Earth-Sun 
libration points. Examples are the Japanese Hiten, the European SMART-1 and 
the NASA GRAIL missions. 

The purpose of this work is to study how the different types of transfer orbits 
to the Moon compare to each other in terms of performance and room occupied in 
the space of the free parameters like the elements of the initial orbit and the launch 
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date. This is done with a Monte Carlo approach, producing a random population 
of transfer orbits within a wide range of initial conditions. 

Section [5] briefly describes the dynamics of the 3- and 4-body problem which 
make the low-energy orbits possible. These are divided into the "outer" type and 
the "inner" type, rather different from each other. Section [3] describes the way the 
numerical simulation was performed, the initial conditions and the algorithm used. 
Section |4] contains an analysis of the Av, the total duration and their dependence 
on the initial conditions of the transfer orbits. Section [5] describes how the data 
could be analysed with the methods of multi-objective optimization to search for 
the preferable type of transfer for a mission. Section [6] states the conclusions. 

2 Types of Earth-Moon Transfers 

2.1 The N-Body Problem 

The study of the dynamics of a small body in the Solar System corresponds to 
searching the solution of a classical, Newtonian, A'' + 1-body problem, where A'' 
is the number of gravitating bodies (Sun and planets) in the system, with given 
initial conditions. 

The problem of flying a probe in the Solar System is that of a point of negligible 
mass that moves along the geodesic defined by its initial conditions in the time 
varying potential given by the revolving planets and the Sun. A precise integration 
of the equations of motion of the probe, even in the simplifying assumption of point 
mass particle, is obtainable by numerical methods once the right hand side (forces) 
is given through high precision ephemerides. The transfer of a spacecraft from the 
Earth to the Moon can be studied, at a high level of accuracy, as that of a body 
of negligible mass in the field of 3 finite mass bodies (Earth, Moon and Sun). 
Analytical solutions are only available in the case of triangular and rectilinear 
configurations. For generic applications, numerical methods are unavoidable. 

The requirements characterizing the PCR3BP are: 

— one of the three bodies is much smaller than the other two and its gravitational 
effect on them is negligible; 

— the two larger bodies, which we will call the primaries A and B travel on 
circular orbits around each other; 

— the third body, which we will call the particle P, moves in the plane of the 
primaries. 

In particular, Joseph Louis Lagrange's work on the 3-body problem led him to the 
discovery of five points of equilibrium, sometimes called libration points, on which 
the whole system's dynamics is based. 

The energy and momentum of the third body (P) are not conserved quantities 
when N = 3. The only first integral known in the PCR3BP is the Jacobi integral 
C, defined as 

C = -2<l>{x,y)-{x^+y^), (1) 

where x and y are the components of the velocity of P and 'P{x, y) is the effective 
potential, containing both gravitational and centrifugal terms. For any given value 
of C, the motion of the particle is spatially confined inside the so called Hill region, 

{{x,y)eR^ : C<-2^{x,y)} , (2) 
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Fig. 1: The a angle and the shape of the gravity gradient induced by the Sun in 
the surroundings of the Earth. The direction of the Sun and the numbers of the 
quadrants (I, II, III, IV) are shown, [(a)] Deviation from the (initially osculating) 
Kepler orbit (dashed gray line) when a is in an even quadrant: the final perigee 
is higher than the initial one. |(b)| When a is in an odd quadrant, the length of 
the major axis is reduced, leading to a lower perigee. 



bounded by zero velocity curves (ZVC), which have different shapes depending on 
the different values of C and of the system's masses. 



2.2 Lunar Transfer and Low Energy Orbits 

In the specific case of a transfer to the Moon, the two primary bodies in the 3-body 
formulation of the problem are the Earth and the Moon, while the spacecraft (P) 
is the third body. This system does not meet the requirements of the PCR3BP 
exactly. However it is close enough to display qualitatively equivalent solutions. 

The traditional method of bringing a payload from the Earth to the Moon is 
through a Hohmann transfer which has a typical total duration Tt of 5 days (from 
LEO). This procedure requires a value of C lower than that of the i*^ and 5*'* 
libration points, thus selecting a configuration of the Hill region where motion is 
allowed in the whole space, without forbidden regions. In this case the duration of 
the flight and the energy of the system allow us to ignore the complex dynamics tied 
to the libration points. Apart from later corrections for the various perturbations, 
this approach only makes use of 2-body celestial mechanics principles. 

Recent research has led to the discovery of another family of traje ctories ex- 
ploiting the properties of t he H i ll region and of t he libration points (iBelbrunol 
Il987l: iBelbruno and Milled. Il990t iKawaeuchi ^Tal Il995t iKoon at all I2OO0L l200lh . 
In order to make this possible, C needs to be close to the critical values, confining 
the space of possible motion and thus allowing the spacecraft to be "captured" 
in the area around the Moon. Some of these trajectories can be achieved only if 
the Jacobi "constant" C varies throughout the flight. The intervention of a fourth 
body is therefore needed to provide an external perturbation to the 3-body system. 
These are sometimes called "low-energy orbits" . 
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Influence of the Sun. The gradient of the solar gravitational field in the surround- 
ings of the Earth is expressed, neglecting the second order terms, with the tidal 
tensor 

^ = ^^[3RR-n\r, (3) 

where R is the Sun-Earth distance vector, u the unit dyadic and r is the Earth- 
spacecraft distance vector. The field introduces a force pushing P away from the 
Earth along the Sun-Earth line and towards it in the direction orthogonal to it, as 
shown in Figure [T] The key parameter relating this effect to the dynamics of an 
ellipt ical orbit is the angl e , q, b etween the Sun-Earth line and the orbit's major 
axis (I Circi and Teofilattol I2OOII) . When a is inside the 2"'^ or 4*'' quadrants in 
the figure, the Sun's gravity field aids the orbit by rising the spacecraft's energy, 
especially when it is close to apogee (Figure[TJi.). The opposite happens when a is 
in one of the other two quadrants, where the spacecraft's energy is lowered (Figure 
Wp.). The overall effect is an increase or decrease of the orbit's major axis growing 
non-linearly with its length. 

It was found ( Belbruno and Milled. I1990I) that it is possible to exploit this 
phenomenon in order to design spacecraft trajectories that reach the Moon with 
lower relative kinetic energy, thus decreasing the overall fuel requirements of the 
mission. The spacecraft is "captured" gravitationally by the lunar field in a process 
called a ballistic capture. 



Ballistic capture. Traditional Hohmann-type transfers approach the Moon with 
a relatively high velocity requiring the use of on-board rockets to slow down the 
motion of P until it becomes bounded to the Moon. In order to obtain a stable, 
bound orbit around the Moon, the energy of the spacecraft with respect to it must 
be decreased to negative values. In the case of the 3-body problem, the condition 
is C > Ci , where C is the Jacobi constant of the spacecraft and Ci is that related 
to the first libration point Li. This corresponds to closing the Li neck of the Hill 
region so that passage between the Earth and Moon domains of possible motion 
is blocked and P is confined in the region around the Moon. 

With ballistic capture, the spacecraft approaches the Moon with a low relative 
velocity and enters directly into a temporarily bound orbit. An example of ballistic 
capture in nature is that of some comets in the combined Sun-Jupiter gravitational 
field. Most famous is comet 39P/Oterma, which was captured for a short time 
around Jupiter in 1936. 

It must be noted that the chaotic 3-body dynamics of the system allows, gen- 
erally, for a temporary (unstable) ballistic capture, although it is thought that 
perr nanent capture ( i.e. the probe remains bound to the Moon forever) can hap- 
pen ( Belbruno . l2004j) . Normally, after the ballistic capture transfer, the spacecraft 
performs one or more revolutions around the Moon before escaping again into 
one of the other realms of possible motion. For this reason an active maneuver is 
necessary to stabilize into a lunar orbit. The clear advantage of using this type 
of trajectory lies in its lower fuel requirements with respect to usual Hohmann 
transfers. 



Outer transfer. The contribution of the Sun described above can be used to natu- 
rally modify the shape of a highly eccentric orbit. This is the case already shown in 
Figure[T] When the apogee distance ra is close to 1.5 x lO^km, the effect becomes 
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strong enough to raise the successive perigee up to about 4 x 10^ km. When this 
happens, the spacecraft can approach the Moon with a lower relative velocity than 
a Hohmann transfer, possibly resulting in a ballistic capture. In this paper this 
type of trajectory is called outer transfer. 

The advantage of this type of transfers, when ballistic capture occurs, is the 
reduced Av in the final phase, when propulsion is needed to stabilize the orbit. 
The data analyzed in the following paragraphs shows that this method results in 
a reduction of up to 44% of this second impulse compared to that of a Hohmann 
transfer. 



Inner resonance transfer. A second type of low-energy transfers does not make 
use of the effect of the Sun. Instead, it is based on the unstable dynamics near the 
Earth-Moon Li libration point. In this case, which we will call inner resonance 
transfer or inner transfer, the apogee distance produced by the initial Av lies 
inside the orbit of the Moon, roughly 3.30 x 10'' km from the Earth. The duration 
of the transfer Tt goes from 85 days upwards, depending on how many Earth orbits 
are completed before falling into the Moon's area of influence. 

In these transfers the effect of the Sun is much weaker and it doesn't help in 
raising the apogee up to the lunar distance, even during the "favorable" a quadrant 
epochs. Yet its contribution can be important when a is in one of the other two 
quadrants, because it can lower the perigee enough to make the spacecraft impact 
the Earth. 

The process by which the capture occurs despite the "short apogee" is the 
execution of one or more resonance hops with respect to the Moon. Each time 
the spacecraft flies by Li its trajectory can shift into a new near-resonant state, 
usually with a longer major axis. This is repeated until the spacecraft is able to 
enter the Moon's sphere of influence. The dynamics behind these resonance hops 
i s tied to the structure o f phase space when P is close in energy and position to Li 
( Belbruno and Marsdenl . Il997t iBelbrunol |2008|) . However its details are still not 
completely understood. 



Interpretation The explana tion given above involving tidal forces and the a angle 
( Circi and Teofilattol l200l[) is only one of several. The original rigorous interpreta - 
tion of the phenomenon is given in IBelbrunol (|l987l) . iBelbruno and Mill^ (|l993l) . 
iBello Mora et all (|2000l) through the concept of weak stability boundary. This is 
the region where the gravitational forces due to the various bodies are similar, 
allowing for ballistic capture and propellant savings. 

Another insightful interpretation involves the stable and unst able manifolds 
related to the periodic orbits at the coUinear libration points (|Conlevl Il968l: 
iKoon et a 1 12OOOI I 200 ih . The phase space position of the spacecraft with respect 
to these two manifolds determines its qualitative behavior in the framework of 
the 3-body problem. By choosing carefully the initial conditions of the orbit it is 
possible to build complex paths going around several celestial objects with few or 
no active maneuvers. 
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3 Method 

The approach followed in this work is statistical, through the search of generic 
transfer orbits from a low Earth parking orbit to a stable lunar orbit. The na- 
ture of the simulation makes it possible to study the complete phenomenology of 
the transfer problem, eliminating any technical or theoretical bias that could be 
included in an ad-hoc algorithm. The aim is to contribute a slightly more com- 
plete picture of the phenomenology of the transfers, rather than discovering new 
methods to create low-energy orbits. 

The simulations are performed using a restricted version of the Solar System 
consisting of Earth, Moon and Sun. They are modeled, in the basic settings, as 
point masses with realistic, non-planar motion, taken from the ephemerides of 
the Jet Propulsion Laboratory's Horizons SysterrQ. The other planets are not 
considered in the present research because they introduce negligible effects in the 
short term. 



3.1 Initial Conditions for the Orbits 

Each transfer is made of two separate ignitions: Avi to begin the transfer at time 
to and Av2, at time t2, to inject the spacecraft into a stable orbit around the 
Moon. The initial condition for all the transfers is a circular low Earth parking 
orbit with a radius ro = 6720 km, and Avi is always applied in a direction parallel 
to the instantaneous velocity 0. 

The ecliptic coordinate system {xec, Dec, Zec) is assumed for the rest of this pa- 
per, unless differently stated. The inclination of the plane of the (circular) parking 
orbit can be varied using two angular parameters: 0, the angle of rotation around 
the and V, the angle of rotation around the yec axis. The angle of rotation 

around the Zec axis is never changed. 

Being the parking orbit circular, the argument of perigee is not defined, and 
the initial conditions for the satellite motion leaving the parking orbit are given 
on the crossing of the parking orbit plane and the ecliptic. This is the only strong 
restriction given to the initial conditions. It greatly reduces the parameter space 
that is explored, but also makes the number of necessary calculations manageable. 

Avi is chosen randomly in an interval between a specified minimum value 
(different choices were tested) and a maximum value that corresponds to the escape 
from the Earth's gravitational well starting from the parking orbit, i.e. Avi = 
3.19 km s~^. After the boost, the trajectory is integrated forward in time for a 
period of 160 days. Among the orbits thus produced, all those that do not enter 
the sphere of radius 10** km around the Moon are discarded. Only those passing 
closer to the Moon than lO"* km without colliding with it are stored. 

After this first screening the algorithm applies a 2- impulse Hohmann maneuver 
to reach a stable 2000 km circular orbit around the Moon (an altitude of 262 km). 
The first impulse of this maneuver is performed at the point closest to the Moon 
of the 160-day orbit. The portion of the trajectory that comes after this point 
of closest approach is discarded. The result is an elliptic orbit with periselenium 

^ http:/ /ssd.jpl.nasa.gov/?horizons 

^ Note that a high trust propulsion capable of instantaneous bursts is assumed. 
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Fig. 2: (a) relative error in energy in the orbital integration for an orbit lasting 
93 days (time is in system's units), ["(b) [ maximum relative error in energy over a 
full orbit, relative to different lengths of the time step. 



distance of 2 x 10^ km. The other impulse is done at periselenium and it circularizes 
the orbit. In total these two maneuvers require a variation of velocity indicated as 

AV2. 



3.2 Orbital Integration 

The integration of the equations of motion of a spacecraft submitted to a variable 
external potential with given initial conditions must be of the highest possible 
precision. Anyway, it is well known that ordinary numerical methods for integrat- 
ing Newtonian equations of motion s beco m e dissipat i ve and exhibit incorrect long 
term behavior (see, e.g., iMenvulJ (|1984 : iMacKavl (|l992l) : ICartwright and Pirc 



(|1992[) ). This is a serious problem performing A'^-body computations, particularly 



when studying the evolution of systems over large time intervals. 

One possible solution is using symplectic integrators, which are a particular 
type of geometric integrators which, respect to ordinary integrators, have the ad- 
vantage to be qualitatively better because they preserve the physical properties 
of the original, although discretized, system's hamiltonian. This allows for a good 
conservation of the characteristics of the dynamical system's time flow properties, 
thing that is particularly appreciated when the simulation is performed over an 
extended time, with an energy error much better constrained than with ordinary 
methods which show an irreversible drift. The advantage of using a symplectic 
algorithm of sufficiently high order is that it gives more reliability of the quality 
of the results at a computer time consumption similar to that of a non-symplectic 
method, so to make it worth using also for integrations limited in time. 

The choice of symplectic methods is wide, for it is possible to construct high 
order integrators ( Yoshida. .1990, ). For the purposes of this paper, we used the 
6*'*-order explicit scheme whose coefficie nts are taken from the first column of the 
Table 1 (SI6A) of iKinoshita et all (|l99l(l . which leads to a conservation of energy 
by a factor 50 better than with the other two possible sets of coefficients reported 
in the same Table. The integration of the spacecraft orbits in the external field was 
done by mean of a properly adapted version of the full N-body code NBSymple, 
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developed bv lCapuzzo-Dolcetta et a The code was also adapted to include 

a multipolar expansion (up to the 4th degree) for the Sun and the Earth potential 
and to account for the influence of the whole set of planets in the solar system, 
but these features were not used for the final computations. 

Although the potential of the Sun is conveniently represented by a truncated 
multipolar expansion, accounting for its deviation from spherical symmetry, it is 
easily seen that the contribution of the quadrupole and hexadecapole (P2 and 
P4 Legendre's terms, respectively) contribute for a fraction of only about 7 x 
10~^^ to the solar gravitational force at one Astronomical Unit. Also the lunar 
potential, at the distances from the Moon of interest for this paper, is sufficiently 
well represented by the monopole term, as can be seen by c omparing i t to hi gher 
moments as deduced from the Lunar Prospector mission bv lHan et aj (|201l|) . So 
in the following calculation we considered, as absolutely sufficient for the scopes of 
this paper, just the solar. Earth and Moon monopole term in their gravitational 
potential. 

The orbit of the point-like spacecraft is integrated with a maximum relative 
error of 2 x 10~^* (see Figure[2]). When the time step is within the range 10"*^ < 
At < 3 X 10~^, our code, working on a standard Intel Xeon CPU is able to generate 
in one minute ~ 250 orbits with a duration of 160 days. 

Further details on the numerical N-bod y code structure and characteristics can 
be found in ICapuzzo-Dolcetta et a il (|2010[) . 



4 Results 

The method employed in the simulation explores a wide portion of the parameter 
space of the Earth-Moon transfer problem. Several kinds of orbits linking the two 
astronomical bodies are found. As it might be expected, most of the trajectories 
found are of limited applicability, with very large Av and long durations. However, 
a significant number of "efi^ective" trajectories, i.e. viable to practical utilization, 
exists. Besides the traditional Hohmann-like transfers, they also include the two 
types of low-energy transfers that exploit the 3- and 4-body dynamics of the Earth- 
Moon-Sun system (inner and outer types, described above). 

A total of 5733 complete transfer trajectories were produced with the following 
4 main simulation modes: 

A. Progression of dates: several transfers each day during the course of one year. 
Lower limit on Avi is such that the minimum ra of the transfer orbit is ra,min = 
3.3 X 10"" km. This includes the Hohmann transfers and the inner and outer 
low-energy transfers. Angles u and are both 0. 

B. Progression of dates: like point 1, but the lower limit on Avi is such that 
ra,min = 8 X 10^ km. This excludes the Hohmann and inner transfers. 

C. progression of values for u\ lower limit on Avi such that ra,min = 8 x 10^ km. 

D. progression of values for 9; lower limit on Avi such that ra,min = 8 x 10"" km. 

The choice of parameters is aimed both at exploring a wide spectrum of possibilities 
and at increasing the number of low-energy orbits to study. Most of the orbits were 
produced with the first two modes. 
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Fig. 3: Left: distribution on the Av\-Av2 plane of 5733 transfers. Green markers: 
simulation mode A (see text); red: mode B; magenta: mode C; blue: mode D. 
Lines corresponding to the equation Avtot = Av\ + Av2 = k for three different 
values of k (as labeled) are shown. Right: relation between ra of an elliptical 
Kepler orbit and Avi, given the initial conditions employed in the simulation. 



4.1 Discussion on the Relation Between Avi and Av2 

A first aspect of the data produced by the simulations is the dependence of the 
second impulse {Av2) on the first (Avi). Avi is only allowed to vary between 
3.08 km and 3.19 km s~^ (the value necessary to reach the escape velocity), 
which results in apogee distances ranging from ra,min = 3.3 x 10'' km to infinity 
(the extreme case being a parabolic orbit). Av2, on the other hand, ranges roughly 
from 0.7 to 2.0 km s"^ 

The left panel of Figure [3] shows a collection of orbits produced for different 
simulation modes. On this plot, an average Hohmann trajectory is located in the 
proximity of the point [3.100 km s~^, 1.093 km s~^]. Different markers denote the 
orbits produced with different modes (see figure caption) . Note that the magenta 
and blue points are not as spread out as the other two modes because they were 
all simulated starting from the same initial date, during a favorable quadrant of 
a. 

The curve in the right panel of Figure |3] shows, for reference, the relation 
between Va and Avi given the initial conditions in use. Note that the Earth-Moon 
distance is 3.85 x 10^ km and the apogee of an outer orbit is ra = 1.2 ~ 1.5 x 10^ 
km. 

Main strip. The distribution in Figure |3] displays a direct positive correlation 
between Avi and Av2 up to Avi = 3.140 km s~^, which corresponds to ra = 7.4 x 
lO'' km. In this region the distribution takes the shape of a strip with increasing 
values of Av2- 

For larger values of Avi , the distribution becomes widely scattered, branching 
also towards lower values of Av2 ■ Most of the transfers on the main strip, excluding 
those in the lower-right area, are inefficient versions of Hohmann transfers and 



10 



R. Capuzzo-Dolcetta, M. Giancotti 



many of them perform several revolutions before coming close to the Moon. The 
reason for this "Hohmann strip" in Figure [3] is straightforward, best explained by 
separating the cases in which Avi is greater than or less than that of the theoretical 
Hohmann which we will here call Avi: 

A. Avi > Avi: here the transfer orbits cross the lunar orbit and go beyond it 
for a while instead of being tangent to it. Trajectories resulting from higher 
Avi 's will therefore approach the Moon with wider angles with respect to the 
lunar instantaneous velocity. The higher this approach angle is, the higher Av2 
needs to be in order to attain a stable final orbit; 

B. Avi < Avi: a traditional Hohmann transfer is not possible because the elliptic 
orbit does not intersect or touch the Moon's trajectory. However the passage 
of P close to the first libration point activates complex 3-body dynamics that 
can lead to a lunar orbit. In particular, most of these orbits perform resonance 
hops, increasing the length of their major axis and eventually falling into a 
ballistic capture around the Moon. The resulting Av2 used to stabilize the 
orbit is usually lower than that of a Hohmann transfer, and in some cases 
similar to that of an outer transfer. 

Distribution for higher Avi. For Avi > 3.140 km s~^, the line of the distribution 
in Figure |3] expands also to low values of Av2- The upper portion is a continuation 
of the previously described Hohmann strip, while the lower branch is composed of 
low-energy outer transfers. 

The difference between the two branches can be entirely ascribed to the 4-body 
dynamics due to the Sun, which becomes important when the point of apogee 
reaches the implied large distances from the Earth. The upper branch is dominant 
when the a angle between the Sun, the Earth and the orbit's apogee is unfavorable, 
i.e. in the l*** or 3'^'' quadrants. Conversely, when a is in the 2"'' or 4*^^ quadrants, 
the probability of obtaining an outer low-energy transfer is higher, producing the 
lower branch. 

All the points of the lower branch around and below the Avtot = 4.0 km s~^ 
diagonal line represent outer low-energy transfers. The lowest value reached is 
Av2 — 0.61 km s~^, that is, 44% less than the Av2 of the average Hohmann 
transfer. 



4.2 Transfer Duration and Av 

Two of the most relevant parameters to be taken into account in the choice of 
the transfer orbit for a mission are the total propellant consumption, expressed 
through Avtot, and the total time of fiight, Tj. Depending on the specific purpose 
of the mission, greater importance may be given to one parameter over the other. 
The propellant used grows exponentially with the Av of a maneuver, and from 
the point of view of the budget it should be kept as low as possible. On the other 
hand, missions that take a long time to reach the target can be unacceptable for 
various reasons, like life support costs in manned missions and related risks. 

The plot in Figure [1] shows the distribution of the orbits in the Avtot-Tt plane. 
The orbits are produced with different choices of the to, v angle and of the lower 
limit to Ta (displayed in the plot with different markers). Figure [3] can be used 
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Fig. 4: 5733 transfers in the Avtot-Tt plane. The same symbols as Figure [3] apply. 



as an aid in the choice of an optimal trajectory given a specific mission, a topic 
explored in more detail in Section [5l 



4.3 Dependence on the Launch Date 

The data obtained from the simulation show that the date of launch influences the 
Av of the transfer orbits. The Moon needs to be in an appropriate part of its orbit 
in order to be reached by the spacecraft, and the required position is different for 
the various types of transfers. The dependence on the lunar cycle is evident in the 
obtained data from the oscillation of the probability of transfer with a period of 
27.3 days, i.e. the period of the Moon's revolution around the Earth. 

The second major effect is given by the angle a between the apogee and the 
Sun (described in Sect. [2}. The gradient of the Sun's gravitational field around the 
Earth's position results in a tidal force acting on the orbiting body, whose direction 
depends on the angular configuration of the Sun-Earth-spacecraft system. 

For short-term orbits, like the traditional Hohmann transfers taking 3 to 5 
days to reach the Moon, the consequences of this force are negligible. However the 
contribution of the Sun cannot be overlooked when considering trajectories with 
ra greater than that of the Moon, with Tt > 80 days. By changing the shape of the 
orbits, the Sun affects the velocity with which the spacecraft passes close to the 
Moon. The result is a variation in the Av2 needed to enter a stable lunar orbit. 

This is visible in Figure [5l where Avtot for different launch dates is shown 
along with the periods in which a is in its various quadrants. Note that the date 
shown for the points is when Avi is applied; a spacecraft in a low-energy orbit 
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Fig. 5: Dependence on the launch date for a set of orbits with 
Avi > 3.14 km s~^. The "unfavorable" quadrants (1^^ and 3'"'^) are shown on a 

darker background. 



comes close to apogee, feeling the effect of the Sun, after roughly 30 to 70 days. 
Therefore, despite the appearance, th e crests and troughs of the distribution do 
correspond to the predicted a phases ( Circi and Teofilattol I2OO1I) . Also, the most 
efficient orbits of the distribution are all located in the favorable quadrants of a, 
as will be shown in Section [5l Figure [TOl 

Of all the transfers with Av2 < 1.0 km s~^, 93% have the apogee angle a 
inside either the 2"^^ or the 4*'' quadrant. The transfers happening when a is 
at the boundary between two quadrants (a = fc7r/2 with k = 1,2,3,4) display 
deformed shapes, because of the different effects of the Sun's gradient along the 
orbit. 

Inner resonance transfers, with ra < 4 X 10^ km do not show a clear dependence 
on a. The mechanism leading to a ballistic capture is, in this case, different and 
mostly independent of the position of the Sun. 



4.4 Low-Energy Orbits 

Existence of low- energy orbits. Taking into account different initial conditions (i.e. 
the A, B, C and D modes described earlier) the number of orbits found to be low- 
energy transfers constitutes a non negligible fraction of the total. The probability 
of producing a low-energy transfer when shooting in a given interval of Avi is 
represented in the histogram in Figure [6l There is no rigorous definition of low- 
energy transfer, but putting an upper limit on the Avtot allows us to exclude 
the Hohmann and the other high-cost transfers. Considering that the average 
Hohmann transfer in this study has Avtot = 4.193 km s~^, we can improvise two 
alternative definitions, one requiring Avtot < 4.0 km s~ and another, stricter one 
requiring Avtot < 3.9 km s~^. 

Figure |6] shows the portion of all the orbits, in each bin of Avi, that fall into 
these two definitions. In both cases the low energy orbits exist in two separate 
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groups, as expected from Section [2l one for low values (inner transfers) and one 
for high values (outer transfers) of Avi. For both definitions there exists a gap 
between the two groups where no conforming transfer takes place. 

The distribution of the transfers in Avi can be translated into ranges of apogee 
distances. The inner transfers thus have apogee distances approximately 3.36 x 
10^ km < Ta < 3.43 x 10^ km. Outer transfers, on the other hand, exist in the 
range 1.23 x 10® km < ra < 1.31 x 10® km. Most of the optimized orbits described 
in other works go further out, up to 1.5 x 10® km from the Earth. The fact that 
these longer orbits are rare in this (non-optimized) simulation is certainly due to 
the greater difficulty in producing them. They occupy a smaller volume in the 
space of the parameters and therefore our Monte Carlo approach is less likely to 
find them. Of course, this more sensitive dependence on the initial conditions is 
the trade-off necessary to obtain higher efficiencies. 

The plots of Figure [7] illustrate the Av distribution of two sets of orbits of 
interest. The "inner set" is composed of all the transfers with Avi < 3.10 km s~^ 
(the value of the mean Hohmann transfer), while the orbits of the "outer set" 
have Avi > 3.140 km s~^. Note that these two sets contain both low-energy and 
high-energy orbits. 

For both sets the majority of the orbits have high Avtot, greater than 4.0 km s~^ 
indicating high-energy orbits. However the left tails of both distributions extend 
down to 3.8 km s~^, where the low-energy orbits are. 

Examples of low-energy orbits. Five different low-energy transfers found in the 
simulation are shown in Figure |8l Orbits [ ^a)| to [(d)] are outer transfers. [ ^a)| has 
i = from the start, while the following three (shown projected on a plane normal 
to the ecliptic) have different initial inclinations. Note that when displayed from 
above the ecliptic, these four transfers all have essentially the same shape. Their 



probability of LET 




Aui (km/s) 



Fig. 6: Histogram showing the portion of orbits that are low-energy, using two 
alternative definitions based on Avtot- 
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outer set, Ar,,,, 





(a) 



(b) 



Fig. 7: Avtot of the inner and outer sets of transfers. The count is normahzed so 
that the sum is unitary, (a) inner set, 499 orbits, Avi < 3.10 km s~^; |(b)|outer 



set, 3538 orbits, Avi > 3.14 km s 



independence to the initial inchnation is due to the nature of the sun's tidal force, 
described in Section [2l constantly pulling towards the ecliptic plane. 

Figures [ ^e)| and [ ^f)| show the same inner transfer in the inertial and in the 
Earth-Moon rotating frame. The sudden changes in the semi-major axis and ec- 
centricity, due to resonance hops, are clearly visible. 



5 Application of Multi-Objective Optimization Methods 

Like most modern day design problems, the choice of an optimal trajectory for 
a space mission has to take into account a number of contrasting objectives, and 
thus belongs to the category called multi-objective optimization (MOO). 

MOO studies the functional relation between the search space, that is the space 
of all the problem parameters, and the objective space, composed by the functions 
that quantify all the object ives. Every objective is translated into one "objective 
function" (Miettineni l200l[ ). 

If there is only one objective, the best solution is the one maximizing or min- 
imizing a suitable objective function, depending on the nature of the problem. 
When the number of objectives is greater than one, however, the conflict between 
them makes some kind of trade-off necessary. The concept of domination of one 
solution over another is introduced to identify which is preferable (or at most 
equal) in terms of all the objectives at the same time. The set of all solutions that 
are not dominated by any other solution, and therefore best unde r at le ast one 
objective, is called Pareto frontier or Pareto-optimal set ( Miettinenl l200lh . 



5.1 MOO for Earth-Moon transfers 

Normally, multi-objective optimization programming starts from the search space 
and uses an algorithm to find the optimal points in the objective space. Such oper- 
ation is not performed in this work because a different approach has been chosen 
for the reasons detailed in the preceding sections. Nevertheless, the production of 
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Fig. 8: (a) Typical outer transfer, i = 0, projection on the x-y plane; [(b)| x-z 
plane projection of another outer transfer, i = 30; (c) outer transfer, i = 90, 
n = 0° , x-z plane; [(d)| outer transfer, i = 90, i7 = 180°, x-z plane; [(e)| |(f)[ inner 
transfer in the inertial and fixed-Earth-Moon frames respectively, i = 28.55, both 

in their x-y planes. 
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Fig. 9: The Pareto front in the Avtot-Tt plane and its 5 order spline 

interpolation. 



a large number of transfer orbits, like the one performed in this work, allows for 
the population of the objective space, where the objective functions are the Avtot 
and the transfer duration Tt- The Pareto set for this problem can be drawn with 
sufficient accuracy, as shown in Figure [9l 

After the identification of the Pareto-optimal transfers, the next step in plan- 
ning the mission is the investigation of the possible best choice of a single solution. 
This task is not easily automated because it depends, obviously, on the specific ob- 
jectives and the restrictions of the mission. Examples of factors to consider would 
be i) the date of launch, ii) whether the mission is manned or unmanned, iii) the 
positions of Sun and Moon at the time of the first impulse, iv) the launch base, v) 
the launcher type and the method of propulsion, etc. All of these may be different 
for each mission project. 

However, the choice of an optimal transfer can be simplified by means of multi- 
objective optimization techniques. In the next subsections this is done for the 2- 
objective problem of minimizing both the transfer duration and the total Av, using 
part of the 7-dimensional search space formed by the initial spatial and velocity 
coordinates and the date of the first impulse. 



5.2 Mapping of the Pareto Front Back to the Search Space 



As a first step we analyse the location of the Pareto-optimal transfers in the space 
of the initial parameters. 
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Fig. 10: Position of the points of the Pareto front (red dots) in the whole 
distribution by date of the orbits. The "unfavorable" quadrants (1^^ and 3'"'^) are 
shown on a darker background. 

If the Pareto front is a continuous curve, it is expected that also its image in 
the search space is a continuous curve. Obviously, the more populated the Pareto 
front is, the more clearly the mapped set will be recognizable as a line. The front 
shown in Figure |9] is determined with 13 points only, and it is obtained from a 
partial exploration of the search space. In these conditions it is anyway possible 
to deduce valuable, although basic, information. 

Figure [10] shows the position of the optimal points (in red) in the Avtot- 
"starting date" plane. The periods of time in which a is in its four quadrants 
is shown in white (2"^^ and 4*^^) and colored (I''* and 3'''^) background. All the 
optimal orbits fall inside the "favorable" (2"*^ and 4*'') quadrants. 

Figure [TT] shows the same plot as Figure [3l that is, the relation between Avi 
and Av2, but with the points corresponding to the Pareto front highlighted in 
red. As detailed in the previous sections, the Avi is directly related to the apogee 
distance. In the plot, the points corresponding to the Pareto set are grouped 
into three distinct areas equivalent to the inner transfers {Avi < 3.1 km i.e. 
apogees inside the orbit of the Moon), the Hohmann transfers {Avi =3.1 km s~^, 
i.e. apogees equal to the Earth-Moon distance) and the outer transfers {Avi > 
3.14 km s~^, apogees longer than 7 x 10^ km). 

The plot in Figure [12] is the distribution, again as a function of the Avi , of 
the distance between the spacecraft and the moon during the capture maneuver 
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A Vi [km/s] 

Fig. 11: The Pareto front (dots joined by tlie red line) in the Avi-Av2 plane. The 
line between the points shows the order in which they appear on the front 
(corresponding to increasing values of Avtot)- The points occupy three distinct 
areas in the plane, highlighted by circles, corresponding to inner, Hohmann and 
outer transfers respectively from left to right. 

Av2- The general distribution of orbits is random, as explained in the previous 
sections of this paper, and all the optimal transfers have capture distances close 
to 2000 km, which is the value that minimizes Av for the final circularization 0. 
This is not surprising, because that the final circularization maneuver contributes 
significantly to the total Av. 

5.3 Scaling Curve 

As was pointed out earlier, there is no algorithm capable of giving a single optimal 
solution for all mission types, and the final decision rests with the humans per- 
forming the optimization. Nevertheless, there are several techniques that aid the 
decision-maker in his choice. One of these is the "weighting method" , in which a 
weight is given to each objective based on its perceived importance. A function ex- 
pressing the distance between the optimal points and a previously selected "ideal" 
or "utopia" point is minimized. Here we use the simple euclidean distance, or L(2) 
norm, between a generic point {xi{p) , X2{p)) and the utopia point {x1,X2)- The 



^ In fact the final orbit around the Moon has a radius of 2000 km. 
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Fig. 12: The Pareto front (red points) in the distribution of distances from the 
Moon at injection, in relation to Av\. 



L(2) norm is defined as: 



(4) 



The Utopia point adopted here is the origin, {Avtot = 0, Tt = 0). Figure [T3l shows 
how different weighting choices deform the front and resuh in different optimal 
points being closest to the utopia point (in normalized units). 

Clearly, the weighting method reduces the problem to a single-objective min- 
imization process. However, it has the drawback of having to choose the weights 
before seeing the respective solution. In some cases a small change in the weights 
can result in large differences in the solution, and the decision-maker has a hard 
time fin ding a clea r answ er. This is partly solved by the "scaling method" devel- 
oped by iKasprzakI ( 200l|) . This method takes each point of the Pareto front and 
computes the weighting (or scaling) that would be necessary to make that solution 
the closest one to the Utopia point. 

The general procedure used in the scaling method is as follows: 

— fit the Pareto optimal set with a high-order polynomial or a spline curve; 

— for each point in the curve, find the local slope m; 

— compute the rescaling necessary to make the slope normal to the distance 
vector that must be minimized; the relation is fc = —m~^. 
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Fig. 13: Examples of three different weiglits for the objective functions and the 
consequent different optimal solutions. The axes are normalized for easier 
handling. Here the xi-axis is deformed for easier viewing, so that the front 
appears to have the same shape, while the circumferences of constant distances 
(in black) appear to be ellipses with different eccentricities. 

The result is a curve relating one of the objective functions to the scaling that 
leads to optimality. For example, by looking at Figure 1131 we see that the curve 
would have (in the normalized units) a value of ki = 0.07 for the red point shown 
in the first graph, k2 = 0.036 for the two points in the second graph and fcs = 0.002 
for the point in the third graph. The full curve, expressed in the original, physical 
units, is shown in Figure [T^ 

The scaling curve obtained in this way facilitates the decision-maker's job, 
because he now has information covering the whole spectrum of solutions, including 
the sensitivity to a change of weighting. The fact that in some cases more than 
one Pareto point corresponds to the same k scaling means that there are multiple 
solutions with the same distance to the utopia point and thus identically optimal, 
like for k2 in Figure [T3l However, the concave part of the front is never optimal in 
the weighting method. Its points never minimize the distance to the utopia point. 
This can be seen in the middle plot in Figure [T31 where the whole arc between 
the two red dots is clearly unreachable by the circumference. All points of this 
kind are shown as a dashed curve in Figure [Tl] and should be excluded from the 
optimization procedure. 

Practical application. The scaling k in Figure [Til can be interpreted physically as 
the number of days of transfer that have, for the specific purposes of the mission 
at hand, the same impact as a Avtot of 1 km/s. 

For example, on a manned mission, 1 day of flight might be as expensive, in 
terms of security and costs, as 1 km/s of the total Av. Using Figure [T4l we look 
for the point in the curve that has fc = 1 days/ (km s~^). Since we have the whole 
scaling curve, our choice of k need not be only one, but can be a range of values. 
In this case we see that the best solution would be for Avtot > 3.91 km/s. From 
this, using the Pareto front, we find that the best transfer orbits are the traditional 
Hohmann transfers. 

On the other hand, transfer duration may have a much lower cost for an un- 
manned mission with science purposes, with a desired scaling ~ 25 days/ (km s~^) 
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Fig. 14: The scaling curve resulting from the Pareto front in Figure [2] and a 
Utopia point at the origin. The dashed part denotes non-optimal solutions. 



In this case the optimal choice would be a low-energy orbit with Avtot — 3.78 km/s 
and (from the Pareto front) Tt ~ 95 days. 

After locating an optimal transfer (or several candidates) on the front, mapping 
back to the search space gives its initial parameters like launch date, inclination 
and Avi. 



6 Conclusions 

In this paper we study the existence and characteristics of low- and high energy 
transfer orbits from the Earth to the Moon in the combined field of Earth, Moon 
and Sun. A wide-ranging and numerically precise study of these orbits contributes 
to a better understanding of the orbital dynamics. Among these transfers, low- 
energy orbits give clear advantages in terms of mission costs. A result of this 
work is that it is easy to obtain successful low-energy orbits just by properly 
setting orbital initial conditions, without resortin g to mid-course maneuvers as, 
for instance, was the case for the Hiten spacecraft ( Belbruno and Milled . Il990l) . 



On the other side, this work has not the scope of a deep understanding of the 
underlying dynamical characteristics of these orbits, which is better studied by 
semi-analytical or geometrical approaches, l ike fo r instance the elegant techniques 
involving invariant manifolds ( Koon et all l200Clt iGomez et aj l2002t iKoon et all 
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l2008h . In particular, it was shown how it is possible using invariant manifol d 
techniques to construct low energy transfer paths to the Moon ( Koon et ail200l|) . 



The existence of such orbits relies strongly on the help of the Sun, and the 
statistics show that a specific range of positions is required for th e Sun for the prac- 
tical realization of low-energy orbits, confirming previous results ( Circi and Teofilattol 



I2OOII) . Another relevant result of this paper is the information about the range of 
apogee distances which most likely produce two-impulse outer low-energy trans- 
fers, which is found to be 1.23 x 10^ km < Va < 1.31 x 10^ km. Moreover, the 
simulations show that the outer orbits are capable of reaching the Moon indepen- 
dently of the inclination of the initial parking orbits, unless they are retrograde. 

As a relevant final point, we discussed the use of multi-objective optimization 
(MOO) methods for the choice of appropriate transfer orbits. The acquired data 
allow us to draw the shape of the Pareto frontier, which is the fundamental element 
of a complete MOO analysis. After determining the Pareto front (set of all the 
optimal points in the objective space) it is possible to map back those points 
onto the search space (that of initial conditions) to get the optimal orbit initial 
conditions (launch date, inclination and first velocity impulse Avi). 

While, clearly, a general technique of resolution of the problem of finding a 
single optimal transfer orbit is not viable because of the variety of contrasting 
objectives, whose specific relevance is highly mission-dependent, we sketched that 
a suitable first approach is that of reducing the MOO to a two-objective problem. 
The two objectives are the total impulse Avtot and the transfer time duration Tt. 
Also in this simplified frame a high degree of indeterminacy remains due to that 
the scope of the mission suggest whether higher relevance is to be given to the 
amount of injected impulse or to the transfer time. 

The final tool used in the selection of a single optimal transfer is the scaling 
curve. The scaling curve shows the best choice given any Tt-Avtot preference fac- 
tor fc, which is mission dependent. We explore two different hypothetical cases 
of manned and unmanned missions. We find, choosing a A: = 1 day/( km s~^) 
for the sake of a manned mission, that optimal orbits are those with Av > 
3.91 km s~^ which coincide with usual Hohmann orbits, while the request of 
fc ~ 25 days/( km s~ ), quite reasonable for unmanned missions, leads to a low 
energy orbit as optimal, with Av ~ 3.78 km/s and Tt ~ 95 days. 
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